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ABSTRACT 

The first gamma-ray burst (GRB) confirmed to be bright enough to be seen 
with the naked eye, GRB 080319B at redshift z = 0.937, allowed for exquisite 
follow-up observations across the electromagnetic spectrum. We present our de- 
tailed optical and infrared observations of the afterglow, consisting of over 5000 
images starting 51 s after the GRB trigger, in concert with our own analysis of 
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the Swift UVOT, BAT, and XRT data. The event is extreme not only in observed 
properties but intrinsically: it was the most luminous event ever recorded at op- 
tical and infrared wavelengths and had an exceedingly high isotropic-equivalent 
energy release in 7-rays. At early times, the afterglow evolution is broadly con- 
sistent with being reverse-shock dominated, but then is subsumed by a forward 
shock at around 1000 s. The overall spectral energy distribution, spanning from 
ultraviolet through near-infrared wavelengths, shows no evidence for a significant 
amount of dust extinction in the host frame. The afterglow evolution, however, is 
highly chromatic: starting at about 1000 s the index shifts blueward before shift- 
ing back to the red at late times. In our deepest late-time observations, we find 
tentative evidence for an optical jet break and a luminous supernova. Finally, 
we examine the detectability of such events with current and future facilities and 
find that such an event could be detected in gamma-rays by BAT out to z — 10.7 
(8a), while the nominal EXIST sensitivity would allow detection to z ~ 32. At 
K band, this source would have been easily detected with meter-class telescopes 
to z ~ 17. 

Subject headings: gamma rays: bursts, gamma-ray bursts: individual: 080319B 



1. Introduction 



The longevity and burst-discovery prowess of the Swift mission (IGehrels et al.l 12004 ) 
has led to a boom in correlative studies of the properties of a large sample of gamma-ray 
bursts (GRBs) observed systematically and uniformly. At the same time, the sheer number of 
bursts encompasses an ever-expanding volume of interesting parameter space, often revealing 
rarities that help to redefine and shape the totality of our understanding of the phenomenon. 
It is in the backdrop of the latter that GRB080319B superlatively reigns. 



At 06:12:49 (UTC is used throughout this paper), GRB 080319B triggered (IRacusin et al. 



2008al ) the Swift Burst Alert Telescope (BAT), the second GRB trigger that day in what 
would be 5 GRB triggers in a 24 hr period. The extreme brightness of the burst at high en- 
ergies, in the X-ray afterglow, and in the ultraviolet /optical/infrared (UV/O/IR) afterglow 
led to a flurry of follow-up observations, many automatically triggered on robotic facilities. 



Contemporaneous imaging, both all-sky and directed, uncovered a fast -rising optica " 
glow which peaked a time from trigger of t gg 18.3 s at V gg 5.3 mag (ICwiok et al 
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2008a 



Covino et al.ll2008l : iSwan et al.ll2008l : iKarpov et al.ll2008al : iSchubel et al.ll2008l : IWozniak et al 
2008 ): this makes the afterglow of GRB080319B the first confirmed counterpart that could 
have been seen with the unaided eye in dark skies. 
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At a redshift of z = 0.937 (jVreeswijk et al.ll2008l ). the event was relatively nearby com- 
pared to the Swift distribu tion of long-duration events (e.g. Jjakobsson et al.ll2006l : lDaigne et al 
20061 : iLe fc Dermerl 120071 ) . yet near the median for redshifts of pre- Swift events. One aim 
of this article is to place the observed properties (particularly energetics) in the context of 
both distributions; we show in §3] that GRB080319B was not only potentially the highest- 
fluence event ever observed, it also had an isotropic-equivalent energy release comparable 
to the highest known values yet recorded. The other significant aim is to analyze in detail 
our long-wavelength observations, in both the optical and IR, with high time cadence on 
moderate-sized robotic telescopes. The quality of the data, coupled with those available in 
the literature and from public archives, reveals a complex evolution of the afterglow that we 
attempt to reconcile with canonical afterglow theory. Unless noted, we assume a concordance 
cosmology with H = 71 km s _1 Mpc -1 , = 0.70, and Q m = 0.3. A redshift of z = 0.937 
corresponds to a luminosity distance of 6011.3 Mpc (distance modulus 43.89 mag). All of 
the results presented herein, though generally consistent with our previous results in GCN 
Circular^], supersede them. 



2. Optical/IR Observations and Reductions 



2.1. PAIRITEL 



The Peters Automated Infrared Imaging Telescope (PAIRITEL; [Bloom et al.l 120061 ). 
responding automatically to the Swift trigg er via an open-source paekaj eonnected to the 
GCN, began taking data on the field at 06:13:39.7 (t — 51 s after the Swift trigger) . The 
1.3 m telescope is equipped with NICMOS3 arrays to simultaneously observe in bands J, 
i7, and K s (1.2, 1.6, and 2.2 /xm, respectively). Each image consists of a 256 x 256 pixel 
array with a scale of 2" pixel -1 . Data are obtained as a double-correlated read, with a 51 ms 
accumulation since reset on the first ("short") read and 7.851 s on the second ("long") read. 
Both the long and the short read exposures are saved to disk. The telescope is dithered 
by ~ V after every third exposure. With no shutter, the standard reductions necessarily 
incor porate per-pixel models o f the combined detector dark current and sky flux fitted over 
time (jWood-Vasey et al.l 120071 ). In a modified version of our automated pipeline, once these 
time-specific frames ("sky+dark") are subtracted from the object frames, they are combined 
into 3—12 file intermediate mosaics with effective integration times of 24—96 s per image, 
which in turn are stacked to form the final mosaics. 



http : 77gcn . gsf c . nasa . gov/ . | 



2 http: //sourcef orge .net/projects/pygcnsock 



-4- 




Fig. 1. — Light curves of the GRB080319B long- wavelength afterglow, fit by our empirical model, 
which allows (and in this case prefers) color change. This is a combin ation of data from the 
GCN Circulars (x symbols, including the prompt light curve as plotted by lKarpov et al.ll2008bl . in 
green), our observations from various ground-based instruments (KAIT, the Lick Nickel lm, and 
PAIRITEL) and our re-reductions of the Swift XJVOT, XRT, and BAT data. The afterglow decays 
extremely rapidly, dropping from mag 5 to 21 in less than one day. For clarity, UV/O/IR data are 
corrected to V [Vega] mag using the model. Individual broken power-law components are shown as 
dotted lines; these are summed together to generate the fitted model (solid line). Different colors 
indicate different filters. Empty points were not used in the fitting in §3.41 
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During the first few minutes on target (253 s in J, 289 s in iif, and 289 s in K s ) ) 
the afterglow was saturated in the exposures, where we have adopted the saturation/non- 
linearity threshold determined from data for the 2MASS South camera^: 37000 counts in a 
single pixel in the J and H bands, and 33000 counts in a single pixel in the K s band. In 
addition, a number of the final science frames suffered from poor sky+dark-frame subtraction 
and could not be photometered or included in the final mosaics. These frames were removed 
following an inspection by eye prior to the construction of the final mosaics. Less than 4% 
of the 7.8 s exposures were removed following this procedure. 



Our initial IR photometry report (IBloom et al.l 120081 ) noted that the afterglow had 



been strongly detected simultaneously in all three filters during the first few minutes of 
observations. Indeed, the ability to detect the transient in single 7.8 s exposures lasted for a 
few hours. In total, we obtained 1822 simultaneous J, H ) and K s observations of the GRB 
(5466 total frames) over the course of 6 hr, though as the afterglow begins to fade we bin 
together individual frames to increase the signal-to-noise ratio and improve readability of 
the light curve. 



2.1.1. Photometry on Unsaturated Frames 

A custom pipeline was used to perform photometry on the 7.8 s exposures constructed 
from the difference of the two readouts which occurred 7.851 and 0.051 s following array reset. 
This pipeline was also used on the intermediate mosaics constructed from the 7.8s frames. 
Zeropoint determinations were made in comparison to the 2MASS catalog. The statistical 
uncertainties are small, while the total uncertainty on any individual image is dominated 
by systematics, especially at early times, associated with our ability to accurately zeropoint 
relative to 2MASS. We quantify this systematic uncertainty empirically by measuring the 
photometric scatter in a mag 10 star (hereafter "SI"), located at a = 14 h 32 m 03 s .55, 5 = 
+36°18 / 29 // .6 (J2000), a few arcmin from the position of the GRB optical afterglow (OA). In 
all three bands the systematic contribution is ^-4%. The resulting light curves are presented 
in Figs. HH2]and listed in Table [U 



3 http : //www . ipac . caltech . edu/2mass/r eleases/allsky/ 
doc/sec4_2 .html . 
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t-ttrigger ( SeC ) 

Fig. 2. — Detailed optical-IR light curve covering much of the first night after the afterglow, 
demonstrating the exquisite time sampling and overall goodness-of-fit of the model. The earliest- 
time PAIRITEL observations are saturated/non-linear and uncertain (see §2.1.11) ; these are not 
included in fitting or modeling. 
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S.i.i?. Photometry on Early- Time Imaging 

Due to the extreme brightness of the early afterglow, a large number of early PAIRI- 
TEL frames are highly saturated or otherwise nonlinear (when the magnitude is < 9). We 
employed two methods to recover the photometry in these early epochs (see Tables EH3]). 
First, since there is sufficient signal in the 51 ms short reads to detect the GRB well above 



the noise when JHK S > 7, we use the procedure described in lEisner et al.l (120071 ) to extract 
the flux from the GRB in the short reads and calibrate this to the 2MASS system using 
stars in the corresponding long reads. The uncertainty of these measurements is primarily 
dominated by the determination of the zeropoint. Second, we developed a simple method 
to extract photometry from saturated "long minus short" frames, using an annulus centered 
on the GRB position. We use SI to determine the zeropoint for each of the individual sat- 
urated long minus short (7.8 s) frames. The inner annulus radius was chosen to exclude 
any saturated pixels in any frame, while the outer annulus radius was selected to provide 
sufficient signal without including too much noise as the wings of the point-spread function 
(PSF) become dominated by noise from the background. The accuracy of these early mea- 
surements of the GRB afterglow is strongly limited by the determination of the zeropoint: 
high near-IR backgrounds and the annulus aperture lead to a large scatter in the determined 
zeropoints for each frame. Thus, the uncertainty for these early measurements is dominated 
by a large systematic term, defined as the scatter of the individual zeropoints. We include 
these points in early-time light-curve plots and in our data table to show the general early 
behavior of the IR afterglow, but do not include them in any fitting or model analysis. We 
emphasize that photometry of these saturated frames is subject to unaccounted-for system- 
atic uncertainties, bu t note qualitative agreement with other submitted early IR imaging 



flRacusin et al.ll2008bh . 



2.2. KAIT and Nickel 



At Lick Observatory, the Katzman A utomatic Imaging Telescope (KAIT: lFilippenko et al. 



200l| ) GRB alert system flLi et al.l 120031 ) responded to GRB080319B automatically. While 
ordinarily the time to slew is < 1 min, since KAIT was following GRB080319A (which trig- 
gered Swift 30 min earlier) our first observations of GRB080319B did not start until 19 min 
after the BAT trigger. We followed the GRB OA with a combination of filters (BVI and 
unfiltered) and with varying exposure times (20 s initially, then 40 s, then 300 s guided expo- 
sures), and continued until 280 min after the BAT trigger. We also observed GRB080319B 
with the Lick 1-m Nickel telescope remotely from the University of California, Berkeley, 
between 116 and 392 min after the BAT trigger, using BVRI filters and exposure times of 
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300 s and 360 s. 



To reduce the KAIT and Nickel data, we used the PSF-fitting technique in IRAF/DAOPHOT. 
Instrumental magnitudes were measured for the GRB OA and several local standard stars, 
and conversion to the standard BVRI system was accomplished using the SDSS calibration 
dCool et alJl2008l : kdelman-McCarthv et allhoosh of the field (the magnitudes in the SDSS 
magnitude system were converted to BVRI following the recipe of Lupton 20050). Two of 
the SDSS standards in our field are shared with the standard star calibration of lHende"n 



(120081 ): because the posted magnitudes for both stars are consistent within the uncertain- 
ties, we adopted the SDSS system due to the much larger number of stars and the ability to 
extend to deep late-time observations. The final KAIT and Nickel photometry is listed in 
Tables @] and El A list of the converted BVRI magnitudes of the bright standards used in 
this calibration is given in Table M 



2.3. UVOT 



To extend the wavelength coverage, we downloaded the Swift UVOT data from the 
quicklook data archive. The Level 2 sky ima ge data in U, B , and V were analyzed according 
to the phot ometry calibration and recipe by lLi et al.l (120061 ). We also followed the procedure 
reported in lPoole et al.l (120081 ). and found that when the GRB was bright, the two procedures 
yield similar results. When the GRB became faint, the Li et al. (2006) procedure yields 
measurements with smaller uncertainties and better overall agreement with the KAIT and 
Nickel data in the B and V passbands, likely due to the smaller adopted photo metry aper 



ture. The Swift UV filters (UVW1, UVM1, and UVW2) were reduced following IPoole et al 
j2008| ). The final UVOT photometry is reported in Table El 



2.4. Gemini Spectroscopy 



An optical spectrum of GRB08031 9B was obtained u nder program GS-2008A-Q-20 



beginning at 08:23 on 19 Mar. 2008 (Cwiok et al.l l2008bh using the Gemini South 8-m 



telescope with GMOS dHook et all 12004 ). We used a slit of width 0.75", the R831 grating, 
and a OG515 filter. Two 1800 s exposures were obtained with slightly different central 
wavelengths of 7000 and 7100 A and read out in 2 x 2 binning. Standard CCD processing and 



4 http: / / www.sdss.org/ dr6 / algorithms / 
sdssUBVRITransform.html#Lupton2005 . 
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spectrum extraction were accomplished with IRAF (for more details, see lFoley et al.ll2006l ). 
The data were extracted using the optimal algorithm of iHornd (119861 ) . The spectrum shows 
a featureless continuum with no strong absorption systems, emission features, or spectral 
breaks. Unfortunately, our wavelength range does not c over the spectral rang e of the VLT 
spectrum where strong absorption features were seen ( Vreeswijk et al.l 2008). With ou r 
spectral coverage, and given the proposed redshift of GRB 080319B (jvreeswijk et al. 2008 ). 
the strongest ISM absorption lines for a galaxy at z — 0.937 are expected to be Ti II 3384 
and the Ca II doublet at A ^ 3950A. We do not detect Ti II A3384 to a 3a rest-frame 
absorption equivalent width limit of W = 0.05 A over a spectral resolution element, while 
the Ca II H & K absorption doublet is blended with the atmospheric A-band absorption. 
No additional features were found in our observed wavelength range of 5950-8150 A. 



2.5. Gemini Imaging 



In the nights following the burst, we began a program of additional imaging using 
Gemini South and Gemini North. On the first night following the event, we acquired 4 x 
180 s of GMOS imaging in each of the g, r, z, and z filters under excellent seeing conditions 
(0.85") despite the low elevation of the target (airmass of 2.5). A second epoch was acquired 
the following night in r only, and additional multicolor epochs were taken later on March 
25 and 28, and April 2. Imaging was taken at Gemini North on April 14 (g, r, i) and 
17 (z). Data wer e reduced using the sta ndard Gemini IRAF package and photometered 
using SExtractor (IBertin fc Arnoutsl Il996) aperture magnitudes. The fie ld was calibrated 
relative to select stars from SDSS DR6 fkdelman-McCarthv et aD I2OO8I . Table [TUP. The 
Gemini filters appear to be somewhat nonstandard compared to the SDSS survey filters, 
so significant color-term corrections were necessary. As the number of comparison stars 
used was large we were able to calculate these individually for each observation. The color 
dependence is about 20% in g on Gemini North and 30% on Gemini South. In r, z, and z 
it is about 10%, 15%, and 5% (respectively) on both telescopes. The final photometry is 
given in Table For use in the light-curve fitting (for which a comparison to the early-time 
BVRI observations is necessary), the gri magnitudes for the afterglow were then converted 
back to VRI using the equations of Lupton (2005). 



2.6. GCN Circulars 



Finally, to supplement our data at very early times and late times, we downloaded 
additional photometry from the GCN Circulars. Of particular note, we downloaded the 
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corrected TORTORA light curve (IKarpov et al.ll2008bl ) containing high-quality optical pho- 
tometry th roughout the prompt phase, and Hu bble Space Telescope (HST) data points at 
late times ( Tanvir et al.ll2008l ; iLevan et al.ll2008l ). which allow us to complete the entire light 
curve of this event. Due to the fact that GCN observations are preliminary and may have 
large calibration offsets, we do not use any of the GCN points in any of our fits or models, 
though we do show them in light-curve plots. Instead, we focus our analysis on the inter- 
mediate phase of the burst, from 300 s to 10 6 s, during which we have good time sampling. 
The GCN data shown in our plots are listed in Table [81 



2.7. Swift BAT and XRT 



Our high-energ y reduction pipeline is desc ribed in detail by iButler et al.l (120071 ) for the 
Swift BAT and bv IButler fc Kocevskil d2007ah for the Swift XRT. GRB080319B exhibits 
one dominant emission episode of duration At « 60 s, composed of multiple unresolved 
spikes. The BAT spectrum in the time interval —1.1 s to 57.4 s is acceptably fit (x 2 /^ = 
13.55/55) by a simple power law with photon index a = — 1.01 ± 0.02 and energy fluence 
(1.96 ± 0.03) x 10 -4 erg cm -2 (15^350 keV). The BAT ca tches only the low -energy portion 



of a spectrum extending beyond £p ea k,obs = 6511^4 keV (IGolenetskii et al.ll2008l ). with an 
isotropic-equivalent energy release in 7-rays of ~ 10 54 erg. There is evidence for a gamma- 
ray tail detected up to ~ 1000 s as seen in Fig. [TJ This extended emiss ion at late times is 
similar to that observe d in the extremely bright BATSE GRB 9 80923 dGiblin et al.lll999h 
and Swift GRB 061007 dMundell et al.ll2007l : Echadv et al.ll2007ah . 



The Swift XR T began obs erving during the tail of the prompt emission phase at t = 66 s. 
We find (see, also. lButlerll2008l ) that the X-ray spectrum is unchanging until quite late times 
(t = 2.91 Ms) despite a break in the X-ray light curve at t ~ 1 ks (Fig. []]). The combined 
XRT/BAT data are well fit by an absorbed power law with T = 1.814 ± 0.011 and an excess 
column density over Galactic of A H = (1.87±0.13) x 10 21 cm -2 at z = 0.937. The early-time 
windowed timing (WT) mode (66 s < t < 4.95 ks) X-ray photon index Ti = 1.814 ± 0.013 
is closely consistent with the late-time photon-counting (PC) mode (4.95 ks < t < 2.91 
Ms) X-ray photon index T 2 = 1.80 ± 0.04, assu ming a constant Nu . From negligible < 1% 
variations in the X-ray hardness ratio (see, e.g., IButler fc Kocevskil l2007bl ) . the magnitude 
of any secular trends in Ti or T 2 must be at the few-percent level or less, and there is no 
evidence for variation in Nr. 
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Fig. 3. — Photometric spectral energy distributions (SEDs) generated over specific time ranges 
during the first night after the burst from combined PAIRITEL, KAIT, Nickel, and UVOT mea- 
surements. All plots have the same horizontal and vertical scale, differing only in the vertical offset. 
The assumed host-galaxy extinction is based on a fit using a combined SED from all epochs (a fit 
assuming no host extinction is also shown for comparison) . A clear red-to-blue transition is evident 
during this time. The three lines show the various fits to the data with the relevant parameters 
noted. "Fixed" indicates that that parameter is held constant in the fit. 
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Fig. 4. — Photometric SEDs generated on subsequent nights from multicolor Gemini South and 
Gemini North observations. The scale and assumed host extinction are the same as in Figure El 
The color is consistent with constant evolution early, but shifts dramatically redward at very late 
times, likely due to the appearance of a very luminous supernova. 




t-Wr ( SeC ) 

Fig. 5. — Evolution of the afterglow spectral index (3 with time. The optical spectral index (3 opt is 
measured from a fit to the UV-optical-IR data from each of 40 different overlapping bins ranging 
from 360 s to 22000 s, plus one additional series of UVOT exposures and the first three multicolor 
Gemini epochs. The X-ray spectral index is also plotted (gray points), with the dashed line showing 
the best-fit value assuming no spectral evolution (which the data are consistent with). The optical- 
to-X-ray spectral index /?ox is defined as the index between the X-ray (normalized at 1 keV) and 
the y-band optical flux. At early times the optical spectral index evolves from red to blue to red 
again. The X-ray, optical, and optical-to- X-ray spectral indices are all consistent at late times, 
which may suggest that both optical and X-ray bands are in the same synchrotron regime. 
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Fig. 6. — Fits to the late-time l^-band data with a combination of an afterglow component 
described with a broken power-law, a constant host galaxy component and a supernova component. 
We assume that the late afterglow, in terms of break-time (11.5 ±3 day) and post-break decay slope 
(2.45 ± 0.35), follows the values derived from the X-ray afterglow (the pre-break decay slope is left 
free to vary, but is strongly constrained by earlier data which is not plotted). The black lines 
represent the sum of all components, directly fit to the data. The blue lines give the afterglow 
plus host galaxy component, whereas the red lines show the supernova contribution. Unbroken 
lines denote the use of the central values (M for Mid) derived from the X-ray fit and a host 
galaxy magnitude of Iq = 26.5. Dotted lines denote the extreme (L)ow case of the (within la 
errors of the X-ray fit) earliest break and steepest post-break decay slope, as well as the faintest 
host galaxy (where we choose conservative errors of 0.5 magnitudes). In this case, at late times, 
the supernova dominates the afterglow. The dashed lines denote the opposite case (H for High), 
a late break and shallow post-break decay as well as a bright host galaxy, which results in the 
smallest supernova contribution, comparable to the afterglow at peak. Even so, we find this "faint" 
supernova component to be 50% more luminous than SN 1998bw at peak. We also plot the three 
z f observations. The two first points show no color evolution, whereas the last point is significantly 
brighter than the extrapolated Iq light curve and gives further evidence of the strong reddening 
due to the rising supernova component. 
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3. Analysis 
3.1. Optical Light Curve 

We fit the optica l data in all filters simultaneously by fitting a series of summed 



1 



Beuermann et al.l ( 119991 ) functions, a generalization of the procedure described by lPerley et al 
(120081 ). Several iterations of different models with varying assumptions and complexity were 
fitted, with our final preferred model motivated by a combination of assumptions of under- 
lying physical behavior (as will be described shortly) and the necessity to fit the data with a 
reasonable x 2 an d without large trends in the residuals. This model contains two elements: 
a monochromatic, rapidly falling component which dominates at early times (< 10 4 s), plus 
a chromatically evolving second component which peaks at about 10 3 - 10 4 s (depending 
on the filter). The first component is actually a sum of three Beuermann functions (this 
was necessary to fit several low-level modulation "wiggles" in the data); the second compo- 
nent is a single Beuermann function, but the peak time is allowed to be a function of the 
filter central wavelength to allow the break (peak) to be chromatic. (Modulations are seen 
at late times as well, but we do not have the temporal coverage to accurately attempt to 
characterize them.) We exclude points after 2 x 10 6 s, which appear to be contaminated by 
additional light (possibly due to a supernova; see §3.7) in most filters. We use only our own 
calibrated data; observations reported in the GCN Circulars are excluded from the fits. Very 
early points from UVOT and PAIRITEL that were heavily affected by pileup, saturation, or 
nonlinearity effects were excluded, and a small number of other conspicuous outliers at late 
times were manually flagged during the fitting process. 

Our model is found to match the data very well, with no obvious residual trends, with 
the exception of a possible chromatic divergence of the very early-time IR measurements from 
the initial power law. The final value for x 2 = 938 (over 592 degrees of freedom) is reasonable, 
if not strictly statistically acceptable, although this is not necessarily a surprise given the 
large number of different instruments and very complicated behavior of the afterglow. A 
modest underestimate of the systematic errors could account for x 2 /dof > 1. All optical 
photometry from our measurements and the GCN Circulars is presented in Figure [TJ Here 
we adopt the convention for the spectral index (/?) and the temporal index (a) such that 
f y oc u~ (3 t~ a . A subset of these observations showing our measurements in more detail is 
shown in Figure El 

The early afterglow decay is extremely rapid: a w 2.24. At about 500 s the light 
cu rve flatte n s sligh tly, and then levels out much more significantly around 1500 s, as noted 



by iLi et al.l (120081 ). Our rapid temporal sampling ends at 20,000 s due to bright morning 



twilight, but we began Gemini observations the following night, and from then until about 
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10 6 s the optical afterglow decays approximately as a simple power law, although some 
limited achromatic variations both above and below the fit may suggest small modulations. 

Our fitting procedure is capable of assigning different colors to different components 
and modeling chromatic breaks, so there is no need to assume achromatic evolution. Indeed, 
attempts at monochromatic fits invariably produced very poor x 2 values and obvious residual 
trends, and color evolution is clear in the time-dependent SED ( §3.21) . Our final preferred 
model is able to fit all the chromatic evolution using the chromatic break of the late-time 
component, which transitions from blue (while rising) to red (while fading). This behavior 
is strongly reminiscent of the predicted evolution of a forward shock at peak, so we fix 
the rising power-law index and the overall color change over the break A(3i^_ a to the 
values predicted by this model. (They are nearly unconstrained by the fits otherwise.) This 
late-time component is allowed to have a different post-break color from the early-time 
component, but the fits suggest they have similar values, with the overall beginning-to- 
end change in spectral index of A/? = 0.006 ± 0.020, despite a notable color change in the 
intermediate region. 



3.2. SED and Extinction Constraints 

Due to the excellent coverage of the photometry across the UV/O/IR spectrum, we can 
strongly constrain the host-galaxy extinction. Our chromatic model can generate SEDs at 
any point in the evolution of the light curve using all available filters, shielding us from the 
possibility that the intrinsic (pre-extinction) SED may not be a simple power law during 
the periods when most of the observations were actually taken (though only to the extent 
to which our model accurately describes the light curve). We can therefore use all available 
filters to fit a single extinction law to the full data in spite of the observed chromatic behavior. 
After correcting for t he small Galactic extinction along this sightline of E(B — V) = 0.01 mag 



flSchlegel et aLlll998h and excluding the Swift UVM2 and UVW2 filters (which are likely to be 
significantly affected by Lyman-a and Lyman-break absorption, respectively), we measure 
host — 0.07 dz 0.06 mag for a fit to the Small Magellanic Cloud (SMC) extinction law. 
Milky Way Galaxy and Large Magellanic Cloud (LMC) laws were also tried, but generally 
had higher values of x 2 — to be sure, the small amount of extinction makes the exact 
choice of extinction law unimportant. Due to the possibility of hydrogen absorption and 
the generally uncertain nature of extinction laws in the far UV, the UVW2 and UVM2 
filters were excluded from these fits. Compared to the X-ray host-galaxy hydrogen column 
of 7V H = (1.87 ±0.13) x 10 21 cm -2 , this suggests a dust-to-gas (actually dust-to-free-metals) 
ratio of about ten times the Galactic value, which is typical for previous GRB sightlines 



(e.g., 
Kann 



Schadv et al 
et al.ll2Q0d ). 



2007bl ) . The small amount of extinction inferred is also typical (e.g., 



We can also generate purely observational time-dependent SEDs. The photometric ob- 
servations taken the night of the burst were divided into 40 overlapping regions, and a simple 
power law was fit to the magnitudes in each region in order to produce contemporaneous 
observed fluxes in each available filter at each region. Galactic extinction, as well as the 
small amount of host extinction correction found above, are removed, and each SED is then 
fit to a power law to estimate the overall UV/O/IR spectral index j3 as a function of time. 
Some representative epochs are shown in Figure [31 These are supplemented by SEDs gener- 
ated directly from the late-time Gemini observations (with a small correction owing to the 
fact that different filters were not observed exactly contemporaneously) , shown in Figure @J 
Finally, we can combine all available data to plot the spectral index f3 as a function of time, 
as shown in Figure El The chromatic behavior of this afterglow is immediately evident. 

3.3. X-ray Fits 

We employ a similar fitting procedure to the X-ray light curve as in the optical and 
distinguish two distinct components: a smooth early decay of a = 1.48 breaks at 2500 s 
to a short lived and poorly constrained (due to lack of observations during this period) 
fast decay of a « 2.8, and then declines as a = 1.40 for the remainder of its evolution 
(possibly except for the last observation in which the afterglow is detected, which appears 
to be significantly below this curve). This differs from the optical value by Aa = 0.17. 
Toward the very end of the XRT observations, the flux appears to decline more rapidly, 
and in our last binned observation is significantly below our fitted model. This appears to 
provide evidence for an additional break at <~ 2 x 10 6 s, but beyond this point the afterglow 
flux faded below the XRT detection threshold and further observations are unavailable, so 
we cannot independently verify this behavior. 

As discussed previously in §2.71 the X-ray spectrum appears invariant over the entire 
burst evolution within fairly stringent constraints. At early times (during the smooth early 
decay) we measure j3x,i — 0.814 ±0.013; during the late phase we measure (3x,2 — 0.80 ±0.04. 

3.4. Testing Canonical Models 

We are able to measure a and (3 very precisely throughout most of the evolution of the 
light curve in both optical and X-ray bands (limited in most cases only by the assumptions 
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of the model itself, e.g., the possible presence of additional underlying components, rather 
than the actual statistical unc ertainty) . In pri nciple this should allow us a strong test of the 
canonical fireball model (e.g. JSari et al.lll998l ). 



During the rapid decay phase of the afterglow, we measure (after extinction correc- 
tion) a opt = 2.2, (3 opt = 0.5 in the UV-optical-IR range (though both values are somewhat 
variable) . Unsurprisingly, this c ombination of lar ge a and low (3 is not well-fit by any canon- 
ical forward-shock model (e.g., IPrice et al.ll2003l ). with the exception of a jet model where 
^opt < v C) which is approximately consistent with both values for p = 2.0 — 2.2. We will 
return to the possibility of an extremely early jet break in §3.61 but the later evolution 
of the light curve (in whi ch it flattens) c asts doubt on this interpretation. Alternatively, 
the reverse-shock model (IKobayashil l2000l ) also predicts a very steep early decay, although 
even this model cannot exactly reconcile the steep decay with the observed spectral index; 
given a = 2.2 this model predicts (3—0.8. Still, this i s a significant improv e ment over any 
forward-shock model. A third possibility, favored by iKumar fc Panaitescul ( 120081 ). is that 
the early steep decay is actually "high latitude" emission: prompt emission whose arrival at 
the observer was delayed because it is slightly off-axis. In this case a value of (3 = a — 2 
= 0.2 is predicted, which is also not consistent with our observations. In fact, the observed 
value is nearly exactly between the reverse-shock and high-latitude values. It is possible 
that both models may contribute in about equal degrees to the observations and produce an 
intermediate spectral index (which may also explain some of the small-scale variation in the 
observed a). 

At about the same time when the optical light curve is rapidly decaying, the X-ray light 
curve is relatively flat, with ax = 1.48. Comparing the relative optical and X-ray fluxes, 
the optical flux overpredicts the X-ray flux at early times; the optical and X-ray bands also 
have different spectral indices, suggesting that a spectral break is present between the two 
bands and that this break is moving upward in frequency with time. Given the disparate 
behavior from the contemporaneous optical light curve and the shift to different behavior at 
late time, the early X-ray emission appears to come from a different origin or process as the 
rest of the light curve; we do not attempt to explain it within the standard model. 

Once the optical light curve begins to level out (in the transition zone between the 
early and late-time components), the spectral index evolves to surprisingly blue values. An 
average SED formed from data points during the extremely well-observed region, using 
observations from K through the UV between 5500 and 12,000 s, constrains the average 
SED to /? opt = 0.07 ± 0.03 (see also Figure [5j). It is impossible to generate (3 < 0.5 during a 
declining light curve within the basic synchrotron model, so this offers strong corroborating 
evidence that the underlying component responsible for the majority of the optical flux is 
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undergoing a break at this time. The canonical model (ISari et al.l Il998l ) could naturally 
explain this portion of the light curve as a characteristic-frequency peak, when v m crosses 
through the UV, optical, and IR bands (in that order), causing a gradual shift in the overall 
color from (3 = — 1/3 to (3 = (p — l)/2. It is in fact difficult to explain this feature in any 
other straightforward way; in a wind med ium the same feature shifts in color only from 



(3 — 1/2 to (3 — p/2 (ILi fc Chevalierl 120031 ) which is never sufficiently blue. Complicated 



prescriptio ns involving numerou s components contributing about equally to the afterglow 



flux (e.g., iRacusin et al.l l2008bl ) may also explain this feature, though at the expense of 



substantial additional complexity. 

Although the optical SED provides good evidence of shifting somewhat back toward the 
red at late times, its overall blueness remains a challenge to explain even after a simple power- 
law decay has set in. This late decay rate is a opt = 1.23±0.02 optically, and ax = 1.30±0.05 
in X-rays. Fitting to a combination of the first three epochs of the late-time Gemini data 
gives (3 opt = 0.51 ± 0.26. If we assume that the X-ray and optical afterglows have the same 
origin, then the X-ray to optical spectral index is /?ox = 0.77 ± 0.01 — consistent with the 
X-ray index itself of (3x = 0.80 ± 0.04. This also forms a strict upper limit on /3 opt ; even if 
a cooling break is present close to the optical band (our observations do not rule out this 
possibility) then (3 opt < /?ox to not underpredict the X-ray flux. 

Given this, and assuming p > 2, we can greatly constrain the available models. The 
jet model (which we favored in the first version of this paper, before late-time data were 
available) is completely ruled out (it requires a > 2.0), as are models where the optical 
band is above the cooling break (which requires (3 > 1.0). This leaves only two possibilities. 
The interstellar medium (ISM) model discussed previously is, in principle, still in agreement 
(barely) with the data: for p = 2.59, a = 1.19 and (3 = 0.79, all within 2a of observed values. 
However this requires that the X-ray regime also be below the cooling frequency {yx < v c ) 
out to extremely late times (at least 10 d), which is not expected physically in the ISM 
model unless the density n and/or the magnetization fraction e# are very low: 

3/2 

< 1. 



cm 



3 J Ux 10- 5 



Alternatively, a wind model with p = 2.0 predicts a opt = 1.25 and (3 opt = 0.50. These are 
reasonable, but the X-ray flux is either greatly over-predicted (if a cooling break is absent, 
we would anticipate /?ox = 0.50 as well), or the X-ray decay is far too fast (if a cooling break 
is present, then ax = 1.0). Either possibility is ruled out to at least 6a. 

We therefore prefer the ISM model to favor all observations. While a wind-stratified 
medium may be made to work with sufficient enhancements (for example, the chromatic 
optical behavior could be explained by an additional emission component, and the late- 
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time discrepancies in the optical and X-rays could be explained by evolution of the cooling 
frequency in a way not predicted by canonical models), this model is simpler as it explains all 
of our observations without the need for such components — with the exception of the very 
early-time data, where standard ass umptions are most expecte d not t o hold. Our conclusio n 
here therefore disagrees with that of iKumar fc Panaitescul (120081 ) and lRacusin et al.l (j2008bl ). 
who prefer a wind-stratified model at late times. 



3.5. Frequency Domain Searches 



Given the brightness of the event, we analyzed the high time cadence 7-ray, X-ray, and 
IR data in search of any significant trends in the frequency do main. The power-density 
spectrum (PDS) for the INTEGRAL data ffeeckmann et alJl2008h . spanning / = 0.03 to 
5 Hz, is well fit by a power law with P(f)df oc f~ 18 . This index here is consistent with the 



—5/3 reported for a number of bright GRBs (IBeloborodov et al J 1 19981 ; IChang fc Yill2000l ). 
and is interpreted as evidence for fully developed turbulence. The PDS from the Swift XRT 
data is essentially (white) noise, apart from 1/f noise at short frequencies. An analysis of 
the first ~2 hr of PAIRITEL J-band photometry also shows a featureless PDS aside from 
1/f noise. 



3.6. Energetics 



The emission from GRB 0803 19B makes it one of the most energetic GRBs ever detected 
and the most energetic thus far seen by Swift . A preliminary analysis of Konus-Wind 



observations reported by lGolenetskii et al.l (120081 ) yields a burst fluence of 5. 72^ 13 x 10 4 erg 



cm 2 in a 20 keV to 7 MeV energy window. This corresponds to a fluence of 5.31 x 10 4 erg 
cm -2 in the BATSE bandpass of 20 keV to 1.8 MeV, making it brighter than the highest 
fluence BATSE burst (4.08 x 10~ 4 erg cm" 2 ; 



equival ent energy release is E\ 



band flAmati et all 12002 ; iBloom et al.ll200lh . 



Kaneko et al.ll2006l ). The rest-frame isotropic 
1.3 x 10 5 4 erg in the standard source frame 1 — 10 4 keV 



If the true collimation-corrected energy (E 7 ) of GRB0803 19B is comparable to the 



median value of E 1 = 1.33 x 10 51 erg found by lBloom et al.l (120031 ). and if we assume average 



values for the efficiency of converting the blast wave 's kinetic energy into gamma rays of 
Ik 



2006 



Kumar et al 



0.5 and the circumburst density n = 1.0 cm ( IGranot fc Kumar 
20071 ). then a jet break is expected at 3.6 x 10 4 s post trigger. A circumburst density greater 



than n — 1 cm 3 only hastens the predicted light-curve break, with n = 100 cm 3 resulting 
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Fig. 7. — Comparison between the observed i?-band light curve of GRB 080319B and those of other GRB 
afterglows, both from t he pr e- Swift as well as the Swift era, shifted to a common redshift of z = 1 with the 
method of iKann et alJ (|2QQ6h . The prompt flash of GRB080319B is clearly shown to be the most luminous 
optical transient ever observed with a high degree of confidence. In spite of this, because of its rapid early 
decay the afterglow at late times is quite unremarkable, and is similar in this regard to the three other 
"ultra-luminous" bursts to date: GRBs 990123, 061007, and 050904. In contrast, the bursts that remain the 
brightest tend to be those with late plateaus and slow decays. 
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in ij et = 7.8 x 10 3 s. Likewise, a collimation-corrected energy less than E 1 = 1.3 x 10 51 erg 
will also push the expected jet break to earlier times, with E 1 = 10 50 erg yielding tj et « 10 3 s. 

No canonical jet break is observed in the GRB light curve. The only unambiguous 
steepening feature is an X-ray break observed at t = 2620 ± 470 s. However, no contempo- 
raneous optical break is observed. There is evidence of a break at very late times (> 10 6 s). 
Unfortunately, the afterglow is extremely faint at this point, but both the last X-ray obser- 
vation, the last g-band detection, and the final HST F606W measurement from the GCN 
circulars are significantly (though only mildly: 2-3a) below our fitted power-law extrapo- 
lation. Furthermore, the fact that the afterglow shifts to dramatically redder colors at this 
point (simultaneously, i and z rise above the fitted curve) suggests that the optical flux may 
be dominated by supernova light (§3.7; host-galaxy light likely contributes as well), and the 
afterglow contribution is actually even smaller than this, which would strongly favor a break. 
Unfortunately, given the lack of available late-time X-ray data we cannot strongly constrain 
whether or not the break is chromatic. 

Treating this feature as a jet break, the very late break time (t = 10 d) imposes strong 
demands on the afterglow energetics. Assuming n = 1 cm -3 and rj 7 = 0.5, this places a 
lower limit on the c ollimation-correcte d energy of 4.0 x 10 51 erg, significantly higher than 



any pre- Swift value (IBloom et al.ll2003l ) and the vast majority of Swift bursts, but compara- 



(Tadiaferri et al. 


|2005 


), 050820A ( 


Cenko et al. 


2006) 



Updike et al 



200 



The lack of an early jet break would then argue that the extreme 
brightness of this GRB is at least partially intrinsic to the explosion itself. 

Alternatively, a jet break may have occurred extremely early (within the first 100 s) 
before the start of our observations. Such a model has been suggested to explain the lack 
of an apparent jet break in anoth er powerful event d etected by Swift , such as GRB 061007 



see 



ISchady et al. 2007a and Mundell et al. 2007 ). In this case, E 1 would actually have 



an unusually low, though not unprecedented, energy release of ~ 10 erg. One could also 
appeal to the minor break in the X-ray light curve (and presume the optical break was hidden 
by other effects), which would indicate a fairly canonical energy near 10 51 erg; either of these 
cases would point to extreme collimation given the large Ej S0 . However, as we have already 
discussed, in these cases the "post-break" light-curve slopes in the opti cal and X-rays o f 
ao = 1.23 and ax = 1.40 are much too shallow for the prediction of t~ p (ISari et al.lll999h . 
where generally we expect the electron index p > 2.0. 

A final possibility is that multiple jets, with v ery different openin g angles, were involved, 
similar to the model proposed for GRB 030329 (IBerger et al.ll2003l ). In this scenario, we 
might expect to see both an early break (from highly collimated, highly relativist ic emission) 
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and a late break (from less collimated, less relativistic emission), presuming that the emission 
at different times was dominated by different jet components. Assuming that the prompt 
emission and optical flash were associated with the narrow jet, the lack of a break at late 
time^l need not imply unusual energetics. The fact that the early-time X-ray break is very 
sharp, and also apparently achromatic, might argue for this interpretation, and the lack of 
a contemporaneous optical break could be explained if the optical counterpart of the early 
X-ray shallo w decay phase were hidden by the bright reverse shock. This model is the one 
favored by iRacusin et al.l (j2008bl ) . 



3.7. Supernova Constraints 



Detailing the nature of the late-time behavior and the evolution of the likely supernova 
counterpa rt is beyond the scope of this paper, but we performed a few basic fits using the 
method of IZeh et al.l (120041 ) of the late-times observations under various model assumptions 
(host brightnesses, using an extrapolation of the late-time light curve and a SN 1998bw-like 
supernova template). In all assumptions we measure a luminosity relative to SN 1998bw in 
the same rest-frame epoch, of greater than one, with an overall distribution of k = 2.3 ± 0.8 
(Figure [6j). In particular, the existence of a late afterglow bre a k and host galaxy with a 
magnitude comparable to the HST measurement of iLevan et al.l (120081 ) favor higher values, 
which provides evidence that this unusually bright GRB may be associ ated with an unusually 
bright (even compared to other GRB supernovae: iFerrero et al.l 120061 ) supernova. For s, the 
relative rest-frame time to maximum light relative to SN 1998bw, we measure s = 0.89±0.10. 



Discussion and Conclusions 



The PDS analysis reflects quantitatively what is qualitatively a light cur ve without 
strong; brightening; featur es. Such "smooth" afterglows have been seen before (e.g.. lLaursen fc Stanek 
20031 : IStanek et al.ll2005l ). yet many GRBs with high-quality optical/IR afterglow observa- 
tions show significant jaggedness, often on timescales less than the time since trigger (e.g., 
Jakobsson et al. 20041 ). The afterglow in this case is not entirely featureless; slow modula- 
tions of less than 10% with At/t are observed at around 200 s and 700 s, but these have 
been seen before in other afterglows, including the undulating afterglow of GRB 030329 



5 To be sure, Racusin et al. ( 2Q08bh claim a late-time jet break but we find no significant evidence of such 
in the XRT data alone. The apparent discrepancy may be due to the different time binnings of the late-time 
data. 
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(ILipkin et all \200i ). 



While an observation of periodicity (/ ^ 1 Hz) in the prompt phase could be considered 
a rea sonable manifestation of a magnetized neutron star powering the initial internal shocks 
(e.g., iMetzger et al.ll2008l ). the absence of such a signature in the INTEGRAL PDS by no 
means rules out such a hypothesis. The concordance of the PDS with other GRBs in the pre- 
Swift sample suggests, prima facie, a similarity between GRB080319B and other events in 
the hydrodynamic properties of the (emitting) outflow; also absent is any substantial evidence 
for a difference in microphysical parameters. We are left with the reasonable conclusion that 
the extreme brightness of GRB080319B has more to do with macroscopic parameters of 
the central engine (in particular, the collimation angle, M e j ect a, initial Lorentz factor, and 
possibly the circumburst medium could all be responsible) than extrema in shock parameters. 

There is a qualitative similarity in the intrinsic behavior of the three brightest afterglow 
events (Fig. [Tj), requiring more than 9 s in the rest frame fro m 7-ray trigger to reach peak 
brightness. Yet, as noted with GRBs 050904 and 990123 feann et all l2007ah . even the 
brightest events fade very rapidly and reach by day 1 a magnitude comparable to that of 
the general population. In our adopted cosmology, M peQj ^ r = —38.3 mag at a time t « 10 
s (rest frame) from trigger. Likewise, M pea ^ r (990123) = —36.0 mag at a time t = 18.1 
s and M peak5r (050904) = -37.1 mag at a time t = 71.2 s. Compared to GRBs 990123 
and 0509 04, GRB080319B does not a ppear to support the proposed brightness vs. i pea k 
relation ( Panaitescu fc Vestrand 20081 ) (i.e., that brighter events are seen to take longer 
from trigger to peak); still, the existence of a rough relation when including fainter events 
suggests that the very brightest events such as these may allow for the longest follow-up 
delay at a given redshift. This afterglow "turn on" delay for some fraction of events has 
important implications for follow-up of high-redshift GRBs since tens of seconds in the rest 
frame becomes a delay of minutes in the observer frame (Fig. [9j). Such a consideration, 
for example, could relax requirements for extremely rapid repointing of satellites to catch 
afterglows at their brightest. 

Figure [8] shows the comparison betw een GRB080319B , quasars, and one of the most 
energetic supernova recorded (SN 2006gv: ISmith et al.l 120071 ). While evolutionary effects in 
all three populations are sure to be important at some level (with QSOs fainter at higher 
redshift, etc.), in the context of probing the high-redshift universe, the overall impression is 
clear: for ~30 min in the rest frame (what would be ~ 4 hr in the observer frame at z = 7), 
GR B 080319b would have been brighter than the brightest known QSO in the universe (see 
also lLazzati et al.ll200ll ; iKann et al.ll2006l for a discussion comparison of the GRB population 
with QSOs). 



We now turn our attention to detectability of such GRBs at high redshift (leaving aside 
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Fig. 8. — Rest-frame comparison of the most luminous optical/IR pro bes of the distant u niverse, 
showing the absolute magnitude (M r , in AB magnitudes as defin ed bvlOke fc Gunnlll982l ) versus 
time of GRB080319B (red circles) and SN2006gy (blu e triangles; 



Smith et al. 



2007^ Transformed 



reen 



light curves of GRB 99 0123 (yellow stars; adapted from lGalama et alJll99d )lmd GRB 050904 (g 
square s; adapted from lKann et al.ll2007bl ) are also shown. For reference, the most luminous known 



QSO (I Schneider et al 



20oJ is shown with a dashed horizontal line; the distribution of SDSS 



QSO magnitudes, adapted from Fig. 6 of I Schneider et all is shown as horizontal banding (darker 
indicates higher density of sources per unit magnitude). The afterglow of GRB080319B was the 
brightest GRB afterglow ever recorded and was at early times ^10 3 times more luminous than the 
most luminous QSO. 
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the question of their existence; see, e.g., iBronim fc Loebl 120071 : iNaoz fc Bromberg 120071 ). 
Given the observed light curve and intrinsic spectrum of this burst (§ 2), to what redshift 
could similar events be detected by present and future missions? The peak photon flux 
from GRB 080319B in 1 s is - 1 x 10 5 9 photons s " 1 . For a BAT threshold flux of 0.8 
photons cm -2 s _1 for i?pk,obs ^ 100 keV ( Bandl 2003). the event would be detected out to 
z = 10.7 (8a). For the nominal EXIST fbrindlavl 120071 ) threshold of 0.2 photons cm" 2 
s _1 , the event would be detected out to z = 32. The ability to detect an afterglow for an 
extreme event of this nature is even more remarkable because of time-dilation effects — 
shifted to z ~ 32 or potentially even further, the afterglow luminosity at late times is nearly 
independent of distance at high redshift. Aside from the effects of line-of-sight absorption 
by neutral hydrodgen, GRB080319B would remain visible even if placed well into the epoch 
of reionization (Fig. [9]). 

We conclude with a rumination on the extrapolation of the features of the afterglow 
toward the low-redshift universe. At z — 0.17, the distance of the nearest non-underluminous 
GRB to date (GRB 030329), this event would peak at R w 1 mag, nearly as bright as the 
brightest stars in the sky. Or, to carry the comparison to its greatest extreme, we might 
envision a situation in which a GRB similar to GRB080319B were to occur in our own 
Galaxy. At a distance of 1 kpc (and neglecting the probable substantial extinction along 
Galactic lines of sight at optical wavelengths), the optical flash would peak at magnitude 
about —28.5, several times the brightness of the Sun! Such an event must assuredly be 
extre mely unusual — the Galac tic GRB rate is probably no greater than 1 per 10 5 to 10 6 



years (IPodsiadlowski et al.l 12004 : see also lStanek et al.ll2006l ). likely only 1% of such bursts 
are collimated toward Earth, and this is among the brightest 0.1% of bursts ever observed. 
Altogether, the rate is probably less than 1 per 10 10 — 10 11 yr: unlikely to have ever happened 
even over the long timescale of geological history, and certainly not a spectacle we can expect 
to witness anytime soon. 
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Fig. 9. — Observability of GRB 080319B in observer-frame K band at a variety of redshifts, based 
on a simple model of the temporal and spectral evolution of the afterglow bootstrapped from the 
y-band light curve. Dotted lines cover the gap in optical (V-band) coverage. The afterglow, which 
would have been as bright as Vega if the event occurred at z — 0.40, remains remarkably bright 
even to z « 16. The purple filled circles show the observed K-band light curve from this paper, 
showing a good agreement with the model of the z — 0.937 (purple) curve (suggesting, too, that the 
interpolation in the V-band time gap is appropriate). The PAIRITEL cumulative 10a point-source 
sensitivity is shown, assuming a nominal start time after GRB trigger of 51 s. For all redshifts 
where the Universe is transparent to Lya photons, this source could have easily been detected by 
sub-meter class telescopes. 
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Table 1. PAIRITEL Observations of GRB 080319B 



^mid 


Filter 


P,vn timp 

l J A I i . L1111C 


Mao- b 


J. 1 l.l-Y 






cop 






344.0 


J 


7.8 


10.389 ± 0.048 


111255.0 ±4815.5 


352.3 


J 


7.8 


10.483 ± 0.048 


102084.5 ± 4419.0 


360.2 


J 


7.8 


10.563 ± 0.048 


94763.2 ± 4102.4 


380.1 


J 


7.8 


10.652 ± 0.048 


87313.1 ± 3779.8 


388.4 


J 


7.8 


10.722 ± 0.048 


81846.4 ± 3544.3 


396.3 


J 


7.8 


10.781 ± 0.048 


77538.9 ± 3357.8 


416.3 


J 


7.8 


10.877 ± 0.048 


70983.9 ± 3074.6 


424.5 


J 


7.8 


10.896 ± 0.048 


69752.5 ± 3021.6 


432.4 


J 


7.8 


10.952 ± 0.048 


66227.7 ± 2868.9 


344.0 


H 


7.8 


9.639 ± 0.042 


142696.2 ± 5421.2 


352.3 


H 


7.8 


9.721 ± 0.042 


132364.7 ± 5029.4 


360.2 


H 


7.8 


9.753 ± 0.042 


128513.3 ± 4883.6 


380.1 


H 


7.8 


9.946 ± 0.042 


107582.1 ±4088.8 


388.4 


H 


7.8 


9.985 ± 0.042 


103750.1 ±3944.2 


396.3 


H 


7.8 


10.115 ±0.042 


92078.9 ± 3501.5 


416.3 


H 


7.8 


10.127 ±0.042 


91066.8 ± 3464.2 


424.5 


H 


7.8 


10.177 ±0.042 


86920.1 ± 3306.4 


432.4 


H 


7.8 


10.245 ± 0.042 


81665.8 ± 3107.6 


344.0 


K s 


7.8 


9.067 ± 0.058 x 


157495.6 ± 8197.4 X 


352.3 


K s 


7.8 


9.308 ± 0.058 x 


126102.6 ± 6564.2 X 


360.2 


K s 


7.8 


9.207 ± 0.058 x 


138421.8 ± 7205.9 X 


380.1 


K s 


7.8 


9.312 ±0.058 


125650.6 ± 6541.1 


388.4 


K s 


7.8 


9.310 ± 0.058 


125949.5 ± 6557.6 


396.3 


K s 


7.8 


9.419 ± 0.058 


113840.4 ± 5927.9 


416.3 


K s 


7.8 


9.508 ± 0.058 


104946.6 ± 5465.3 


424.5 


K s 


7.8 


9.546 ± 0.058 


101331.5 ± 5277.5 


432.4 


K s 


7.8 


9.638 ± 0.058 


93031.2 ± 4845.6 



Note. — Because of the very large number of exposures acquired, only a few representative points 
are given. The full table of photometry containing all 406 points is available online. Includes only 
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non-saturated exposures; photometry from the saturated epochs is given in Tables [2] and El 
a Exposure mid-time, measured from the Swift trigger (UTC 06:12:49). 
b Observed value; not corrected for Galactic extinction. 
x Point not used in modeling. 
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Table 2. PAIRITEL Short-Read Observations of GRB 080319B During Saturated Regime 



^mid 

sec 


Filter 


Exp. time 

sec 


Mag. b 




Flux b 


59.3 


J 


0.051 


5.206 ± 0.055 


1.32 x 


10 7 ± 6.52 x 10 5 


67.3 


J 


0.051 


5.627 ±0.055 


8.94 x 


10 6 ± 4.43 x 10 5 


87.3 


J 


0.051 


6.470 ± 0.055 


4.11 x 


10 6 ± 2.04 x 10 5 


95.5 


J 


0.051 


6.847 ± 0.056 


2.91 x 


10 6 ± 1.45 x 10 5 


103.5 


J 


0.051 


6.954 ± 0.056 


2.63 x 


10 6 ± 1.32 x 10 5 


59.3 


H 


0.051 


4.520 ± 0.052 


1.59 x 


10 7 ± 7.47 x 10 5 


67.3 


H 


0.051 


4.853 ± 0.052 


1.17 x 


10 7 ± 5.50 x 10 5 


87.3 


H 


0.051 


5.713 ± 0.052 


5.30 x 


10 6 ± 2.50 x 10 5 


95.5 


H 


0.051 


6.101 ±0.052 


3.71 x 


10 6 ± 1.75 x 10 5 


103.5 


H 


0.051 


6.225 ± 0.052 


3.31 x 


10 6 ± 1.56 x 10 5 


67.3 


K s 


0.051 


4.362 ± 0.048 


1.20 x 


10 7 ± 5.20 x 10 5 


87.3 


K s 


0.051 


5.192 ± 0.048 


5.59 x 


10 6 ± 2.43 x 10 5 


95.5 


K s 


0.051 


5.308 ± 0.048 


5.02 x 


10 6 ±2.18 x 10 5 


103.5 


K s 


0.051 


5.675 ± 0.048 


3.58 x 


10 6 ± 1.56 x 10 5 



a Exposure mid-time, measured from the Swift trigger (UT 06:12:49). 
b Observed value; not corrected for Galactic extinction. 
x No saturated values were used in modeling. 



Note. — Because of the large number of exposures acquired, only a few representative points are 
given. The full table of photometry containing all 406 points is available online. 
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Table 3. PAIRITEL Annulus Photometry of GRB 080319B During Saturated Regime 



^mid 


Filter 


Exp. time 


Mag. b 




Flux b 


sec 




sec 








55.0 


J 


7.8 


4.595 ± 0.133 


2.31 x 


10 7 ± 2.66 x 10 6 


63.2 


J 


7.8 


5.070 ± 0.133 


1.49 x 


10 7 ± 1.72 x 10 6 


71.2 


J 


7.8 


5.560 ±0.133 


9.51 x 


10 6 ± 1.10 x 10 6 


91.2 


J 


7.8 


6.572 ±0.133 


3.74 x 


10 6 ± 4.31 x 10 5 


99.4 


J 


7.8 


6.710 ±0.133 


3.30 x 


10 6 ± 3.80 x 10 5 


107.4 


J 


7.8 


7.004 ± 0.133 


2.51 x 


10 6 ± 2.90 x 10 5 


55.0 


H 


7.8 


4.387 ±0.191 


1.80 x 


10 7 ± 2.90 x 10 6 


63.2 


H 


7.8 


4.805 ±0.191 


1.22 x 


10 7 ± 1.97 x 10 6 


71.2 


H 


7.8 


5.159 ±0.191 


8.84 x 


10 6 ± 1.42 x 10 6 


91.2 


H 


7.8 


5.818 ±0.191 


4.82 x 


10 6 ± 7.76 x 10 5 


99.4 


H 


7.8 


6.386 ±0.191 


2.86 x 


10 6 ± 4.60 x 10 5 


107.4 


H 


7.8 


6.421 ±0.191 


2.76 x 


10 6 ± 4.45 x 10 5 


55.0 


K s 


7.8 


3.853 ± 0.463 


1.92 x 


10 7 ± 6.65 x 10 6 


63.2 


K s 


7.8 


4.200 ± 0.463 


1.39 x 


10 7 ± 4.83 x 10 6 


71.2 


K s 


7.8 


4.510 ±0.463 


1.05 x 


10 7 ± 3.63 x 10 6 


91.2 


K s 


7.8 


5.130 ± 0.463 


5.91 x 


10 6 ± 2.05 x 10 6 


99.4 


K s 


7.8 


5.659 ± 0.463 


3.63 x 


10 6 ± 1.26 x 10 6 


107.4 


K s 


7.8 


5.812 ±0.463 


3.16 x 


10 6 ± 1.09 x 10 6 



Exposure mid-time, measured from the Swift trigger (UT 06:12:49). 
'Observed value; not corrected for Galactic extinction. 
No saturated values were used in modeling. 



Note. — Because of the large number of exposures acquired, only a few representative points are 
given. The full table of photometry containing all 406 points is available online. 



-37- 



Table 4. KAIT Observations of GRB 080319B 



*mid a Filter Exp. time 


Mag. b 


Flux b 


sec sec 







1217.0 


clear 


20.0 


14.116 ±0.008 


7008 


1 


±51.4 


1308.0 


clear 


20.0 


14.261 ± 0.007 


6132 





±39.4 


1404.0 


clear 


20.0 


14.424 ± 0.011 


5277 


2 


±53.2 


1495.0 


clear 


20.0 


14.548 ± 0.013 


4707 


6 


±56.0 


1586.0 


clear 


20.0 


14.651 ±0.011 


4281 


5 


±43.2 


2098.0 


B 


20.0 


15.623 ± 0.045 


2327 





±94.5 


2249.0 


B 


20.0 


15.750 ± 0.041 


2070 


1 


±76.7 


2373.0 


B 


20.0 


15.814 ± 0.042 


1951 


6 


±74.1 


2546.5 


B 


40.0 


15.911 ±0.037 


1784 


8 


±59.8 


2791.0 


B 


40.0 


16.094 ± 0.035 


1508 





±47.8 


1157.0 


V 


20.0 


14.345 ± 0.022 


6745 


3 


± 135.3 


1248.0 


V 


20.0 


14.442 ± 0.022 


6168 


8 


± 123.7 


1339.0 


V 


20.0 


14.651 ± 0.022 


5088 


6 


± 102.1 


1435.0 


V 


20.0 


14.815 ± 0.025 


4375 


2 


±99.6 


1526.0 


V 


20.0 


14.939 ± 0.026 


3903 





±92.4 


1188.0 




20.0 


13.777 ± 0.013 


7502 


4 


±89.3 


1279.0 




20.0 


13.879 ±0.012 


6829 


7 


±75.1 


1370.0 




20.0 


14.016 ±0.015 


6020 





±82.6 


1465.0 




20.0 


14.176 ±0.020 


5195 


2 


±94.8 


1557.0 




20.0 


14.328 ± 0.018 


4516 


5 


±74.3 



Note. — Because of the very large number of exposures acquired, only the first five exposures in 
each band are given. The full table of photometry is available online. 

a Exposure mid-time, measured from the Swift trigger (UTC 06:12:49). 

b Observed value; not corrected for Galactic extinction. 
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Table 5. Nickel Observations of GRB 0803 19B 



^mid 


Filter 


Exp. time 


Mag. b 


Flux b 


sec 




sec 






7154.7 


B 


360.0 


-1 t7 A /"» T" 1 A AAr 

17.365 ± 0.005 


467.7 ± 2.1 


8579.2 


B 


360.0 


17.646 ± 0.006 


361.1 ± 2.0 


-i a a -i a a 

10013.4 


B 


360.0 


17.866 ± 0.008 


AA /I A 1 A A 

294.8 ± 2.2 


11512.3 


B 


360.0 


i c\ r\r\r\ i a aaa 

18.090 ± 0.008 


A A A A 1 -1 A 

239.9 ± 1.8 


-i a a -i f a 

13015.9 


B 


360.0 


-i a a a a i a a -i -i 

18.248 ± 0.011 


A Af7 A 1 A -1 

207.4 ± 2.1 


14473.6 


B 


360.0 


-i a a a a i A a -i -i 

18.389 ± 0.011 


-| A A "I 1 "1 A 

182.1 ± 1.8 


16436.8 


B 


360.0 


-i a rnf7 i a a -i k 

18.587 ± 0.015 


151.8 ± 2.1 


17937.4 


B 


360.0 


-1 A l~7l—7f\ 1 A A "1 A 

18.770 ± 0.018 


-1 A A A 1 A -1 

128.2 ± 2.1 


-i /"v a a a i 

19398.1 


B 


360.0 


-i o rife i a a -i r" 

18.858 ± 0.015 


-i -i a a i -1 r* 

118.2 ± 1.6 


20965.3 


B 


360.0 


18.970 ± 0.023 


-i r\ r* 1-7 \ A A 

106.7 ± 2.2 


1-7 A A a a 

7444.0 


V 


300.0 


~\ P7 C\T" ~\ 1 A AAf 

17.251 ± 0.005 


464.1 ± 2.1 


8860.5 


V 


300.0 


~\ <~7 A O A 1 A f\f\<~7 

17.484 ± 0.007 


O ^7 /I f 1 A | 

374.5 ± 2.4 


10291.9 


V 


300.0 


-| 1-7 r7A^7 1 A a r\ i—r 

17.727 ± 0.007 


AAA A 1 "1 A 

299.4 ± 1.9 


H829.3 


V 


300.0 


-I OA A 1 A r\C\>~7 

17.899 ± 0.007 


255.5 ± 1.6 


13294.2 


V 


300.0 


-1 A A/1 A 1 A AAf7 

18.062 ± 0.007 


A -1 A A 1 -1 A 

219.9 ± 1.4 


14767.2 


V 


300.0 


-1 A A A /» 1 A A1 /I 

18.236 ± 0.014 


1 A 1-7 A 1 A A 

187.3 ± 2.4 


16714.4 


V 


300.0 


-io in /I i a a -i a 

18.424 ± 0.012 


157.5 ± 1.7 


"1 A A A /"* f 

18236.5 


V 


300.0 


-i a r— i— r\ i a a -i r* 

18.559 ± 0.016 


-1 A A -1 1 A A 

139.1 ± 2.0 


19700.4 


V 


300.0 


-i a /iof7 i a a -i a 

18.637 ± 0.013 


-1 AA t" 1 -1 

129.5 ± 1.5 


aaaa/"* a 

22996.8 


V 


300.0 


l o AAA 1 A AA"1 

18.902 ± 0.021 


-I A -1 /I 1 -1 A 

101.4 ± 1.9 


7782.0 


R 


300.0 


16.991 ± 0.006 


496.1 ±2.7 


9196.4 


R 


300.0 


17.239 ± 0.005 


394.8 ± 1.8 


10635.4 


R 


300.0 


17.435 ± 0.008 


329.6 ± 2.4 


12181.0 


R 


300.0 


17.637 ± 0.009 


273.7 ±2.3 


13652.0 


R 


300.0 


17.799 ± 0.007 


235.7 ± 1.5 


15110.2 


R 


300.0 


17.917 ±0.007 


211.4 ±1.4 


17071.9 


R 


300.0 


18.107 ±0.011 


177.5 ± 1.8 


18582.1 


R 


300.0 


18.218 ± 0.008 


160.3 ±1.2 


20035.8 


R 


300.0 


18.308 ± 0.017 


147.5 ± 2.3 


23331.8 


R 


300.0 


18.541 ± 0.021 


119.0 ±2.3 


8122.0 


I 


300.0 


16.797 ± 0.004 


464.7 ± 1.7 
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Table 5 — Continued 



t mid a Filter 


Exp. time 


Mag. b 


Flux b 


sec 


sec 




(iJy 


9531.2 I 


300.0 


17.013 ± 0.005 


380.9 ± 1.7 


10995.4 I 


300.0 


17.189 ±0.005 


323.9 ± 1.5 


12520.7 I 


300.0 


17.394 ± 0.005 


268.2 ±1.2 


14003.3 I 


300.0 


17.537 ± 0.009 


235.1 ± 1.9 


15973.2 I 


300.0 


17.704 ± 0.008 


201.6 ±1.5 


17412.9 I 


300.0 


17.817 ±0.011 


181.6 ± 1.8 


18936.7 I 


300.0 


17.921 ±0.010 


165.0 ±1.5 


20447.9 I 


300.0 


18.045 ± 0.014 


147.2 ± 1.9 


23665.5 I 


300.0 


18.111 ±0.118 


138.5 ± 14.3 



a Exposure mid-time, measured from the Swift trigger (UTC 06:12:49). 
b Observed value; not corrected for Galactic extinction. 
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Table 6. UVOT Observations of GRB 080319B 



t • a 

^mid 


Filter 


Exp. time 


Mae: b 


Flux b 


sec 




sec 






6156.0 


UVW1 


196.6 


16.146 ±0.043 


343.2 ± 13.3 


17398.0 


UVW1 


885.6 


17.645 ± 0.042 


86.3 ±3.3 


28047.0 


UVW1 


427.9 


18.316 ± 0.086 


46.5 ±3.5 


45393.0 


UVW1 


427.9 


18.909 ±0.121 


26.9 ±2.8 


75628.0 


UVW1 


200.6 


19.266 ± 0.213 


19.4 ±3.5 


5542.0 


UVW2 


196.6 


16.192 ±0.046 


239.6 ± 9.9 


10709.0 


UVW2 


885.6 


17.231 ± 0.036 


92.0 ± 3.0 


22273.0 


UVW2 


266.3 


18.394 ±0.116 


31.5 ±3.2 


28931.0 


UVW2 


856.1 


18.631 ± 0.073 


25.3 ±1.6 


46279.0 


UVW2 


856.1 


19.228 ±0.102 


14.6 ±1.3 


76050.0 


UVW2 


401.3 


19.888 ± 0.220 


8.0 ±1.5 


5951.0 


UVM2 


196.6 


15.923 ± 0.052 


284.4 ± 13.3 


12528.0 


UVM2 


334.6 


17.144 ±0.071 


92.4 ±5.8 


16491.0 


UVM2 


885.6 


17.493 ± 0.052 


67.0 ±3.1 


30028.0 


UVM2 


611.2 


18.381 ±0.101 


29.6 ±2.6 


47376.0 


UVM2 


610.0 


19.039 ±0.144 


16.1 ±2.0 


631.0 


U 


19.4 


12.474 ± 0.028 x 


19697.0 ± 501.5 X 


785.0 


U 


19.4 


12.968 ± 0.030 


12496.8 ± 340.6 


1430.0 


U 


19.4 


14.108 ±0.042 


4373.2 ± 165.9 


1591.0 


U 


19.5 


14.342 ± 0.046 


3525.3 ± 146.2 


4926.0 


u 


196.6 


16.015 ± 0.029 


755.1 ± 19.9 


6362.0 


u 


196.6 


16.329 ± 0.033 


565.5 ±16.9 


18304.0 


u 


295.0 


17.864 ± 0.056 


137.5 ±6.9 


18609.0 


u 


231.2 


17.886 ± 0.064 


134.8 ± 7.7 


28486.0 


u 


213.4 


18.498 ± 0.089 


76.7 ±6.0 


45833.0 


u 


213.4 


19.435 ±0.148 


32.4 ±4.1 


656.0 


B 


9.6 


13.334 ± 0.042 x 


19160.2 ± 727.0 X 


811.0 


B 


9.6 


13.913 ± 0.046 


11240.9 ± 466.3 


1454.0 


B 


19.4 


15.012 ± 0.046 


4085.1 ± 169.5 


1616.0 


B 


19.5 


15.291 ± 0.052 


3159.4 ± 147.7 


5131.0 


B 


196.6 


16.962 ± 0.033 


678.0 ± 20.3 


6567.0 


B 


196.6 


17.350 ± 0.039 


474.2 ± 16.7 
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Table 6 — Continued 





Filter 


Exp. time 


Mag. b 


Flux b 


sec 




sec 




(iJy 


28708.0 


B 


213.4 


19.417 ±0.112 


70.7 ±6.9 


46055.0 


B 


213.3 


20.036 ±0.171 


40.0 ± 5.8 


75943.0 


B 


100.2 


20.441 ± 0.316 


27.5 ±6.9 


711.0 


V 


19.5 


13.373 ± 0.040 x 


16512.0 ± 597.2 X 


975.0 


V 


393.5 


14.337 ± 0.013 x 


6795.2 ± 80.9 X 


1518.0 


V 


19.4 


15.025 ± 0.076 


3605.8 ± 243.8 


1680.0 


V 


19.5 


15.127 ±0.081 


3282.5 ± 236.0 


5746.0 


V 


196.6 


16.955 ± 0.058 


609.5 ± 31.7 


11616.0 


V 


295.0 


17.948 ± 0.078 


244.2 ± 16.9 


11919.0 


V 


295.1 


17.895 ± 0.078 


256.4 ± 17.8 


12222.0 


V 


295.0 


18.041 ± 0.082 


224.2 ± 16.3 


29806.0 


V 


213.4 


19.491 ± 0.256 


59.0 ± 12.4 



a Exposure mid-time, measured from the Swift trigger (UTC 06:12:49). 
b Observed value; not corrected for Galactic extinction. 
x Point not used in modeling. 
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Table 7. Gemini Observations of GRB 080319B 



^mid 


Filter 


Exp. time 


Mag. b 


Flux b 


sec 






sec 




(iJy 


89309 


g 


5 


X 


180 


20.670 ± 0.100 


19.59 ± 1.72 


90356 


r 


5 


X 


180 


20.520 ± 0.060 


22.49 ± 1.21 


91307 


i 


5 


X 


180 


on oon i r\ r\r r\ 

20.380 ± 0.050 


or rn i i i r 

25.59 ± 1.15 


92470 


z 


5 


X 


180 


20.310 ± 0.050 


27.93 ± 1.26 


174960 


r 


5 


X 


180 


21.510 ±0.060 


9.037 ± 0.486 


520390 


g 


6 


X 


180 


23.450 ± 0.190 


1.514 ±0.243 


521589 


r 


6 


X 


180 


23.340 ± 0.090 


1.675 ±0.133 


522792 


i 


6 


X 


180 


22.900 ± 0.060 


2.512 ±0.135 


523990 


z 


6 


X 


180 


22.940 ± 0.090 


2.477 ±0.197 


780120 


r 


6 


X 


180 


23.670 ± 0.070 


1.236 ± 0.077 


781560 


i 


6 


X 


180 


23.280 ± 0.060 


1.770 ± 0.095 


1209960 


g 


6 


X 


180 


24.440 ± 0.080 


0.6081 ± 0.0432 


1211220 


r 


6 


X 


180 


24.100 ± 0.080 


0.8318 ± 0.0591 


1212480 


i 


6 


X 


180 


23.810 ± 0.070 


1.086 ± 0.068 


2286720 


g 


6 


X 


180 


25.860 ±0.110 x 


0.1644 ±0.0158 x 


2287800 


r 


6 


X 


180 


25.050 ± 0.090 x 


0.3467 ± 0.0276 x 


2285280 


i 


6 


X 


180 


24.360 ± 0.080 x 


0.6546 ± 0.0465 x 


2531160 


z 


3 


X 


180 


23.940 ± 0.130 x 


0.9863 ±0.1113 x 



a Exposure mid-time, measured from the Swift trigger (UTC 06:12:49). 
b Observed value; not corrected for Galactic extinction. 
x Point not used in modeling. 
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Table 8. GCN Observations of GRB 080319B 



sec 


Filter 


Exp. time 

sec 


Mag. b 




Flux b 






-11.0 


clear 


10.0 


> 11.48 


< 7.94 x 10 4 


Pi c 


GCN 7445 d 


3.0 


clear 


10.0 


9.83 


3.63 x 


10 5 


Pi 


GCN 7445 


17.0 


clear 


10.0 


5.76 


1.54 x 


10 7 


Pi 


GCN 7445 


32.0 


clear 


10.0 


6.00 


1.24 x 


10 7 


Pi 


GCN 7445 


79.0 


clear 


10.0 


8.26 


1.54 x 


10 6 


Pi 


GCN 7445 


93.0 


clear 


10.0 


8.77 


9.64 x 


10 5 


Pi 


GCN 7445 


108.0 


clear 


10.0 


9.10 


7.11 x 


10 5 


Pi 


GCN 7445 


122.0 


clear 


10.0 


10.27 


2.42 x 


10 5 


Pi 


GCN 7445 


163.0 


clear 


10.0 


10.50 


1.96 x 


10 5 


Pi 


GCN 7445 


177.0 


clear 


10.0 


11.10 


1.13 x 


10 5 


Pi 


GCN 7445 


252.0 


clear 


10.0 


11.21 


1.02 x 


10 5 


Pi 


GCN 7445 


296.0 


clear 


10.0 


11.79 


5.97 x 


10 4 


Pi 


GCN 7445 


310.0 


clear 


10.0 


11.95 


5.15 x 


10 4 


Pi 


GCN 7445 


37.2 


clear 


10.0 


5.60 


1.79 x 


10 7 


RAPTOR 


GCN 7464 e 


24.6 


clear 


5.0 


5.35 


9 9^ v 


in 7 


ROTSE 


GCN 7470 f 


111.7 


clear 


20.0 


8.49 


1.25 x 


10 6 


ROTSE 


GCN 7470 


8.9 


clear v 


- 


8.198 ±0.579 


1.94 x 


10 6 ± 8.02 x 10 5 


TORTORA 


GCN 7558§ 


10.1 


clear v 


- 


6.967 ±0.188 


6.03 x 


10 6 ± 9.58 x 10 5 


TORTORA 


GCN 7558 


11.5 


clear v 


- 


6.650 ±0.141 


8.07 x 


10 6 ± 9.83 x 10 5 


TORTORA 


GCN 7558 


60300.0 


R 




19.80 


3.73 x 


10 1 


Canarias 


GCN 7476 h 


55830.0 


clear 


1500.0 


19.080 ± 0.020 


7.24 x 


10 1 ± 1.32 x 10° 


AAVSO 


GCN 7484 1 


70500.0 


clear 


2040.0 


19.000 ±0.100 


7.80 x 


10 1 ± 6.86 x 10° 


AAVSO 


GCN 7484 


74831.0 


R 




20.000 ± 0.300 


3.10 x 


10 1 ± 7.50 x 10° 


Brno 


GCN 7504J 


43934.4 


R 




19.10 


7.11 x 


10 1 


Kharkiv 


GCN 7519 k 


46270.8 


R 




19.90 


3.40 x 


10 1 


Kharkiv 


GCN 7519 


71208.0 


R 




20.10 


2.83 x 


10 1 


Kharkiv 


GCN 7519 


1656000 


R n 




24.35 


5.65 x 


10" 1 


HST 


GCN 7569 1 


1656000 


1° 




23.76 


7.62 x 


10" 1 


HST 


GCN 7569 


4590000 


R n 




26.280 ±0.100 


9.55 x 


10" 2 ± 8.40 x 10" 3 


HST 


GCN 7710 m 


4590000 


1° 




25.460 ±0.100 


1.59 x 


10" 1 ± 1.40 x 10" 2 


HST 


GCN 7710 
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Note. — Most of the TORTORA measurements have been omitted to save space. The full table 
of GCN photometry used in our plot is available online. 

a Exposure mid-time, measured from the Swift trigger (UTC 06:12:49). 

b Observed value; not corrected for Galactic extinction. 

c "Pi-of-the-Sky" Collaboration 



^Cwiok et al.1 fj2008ah 



Wozniak et all (120081 ) 



Swan et all (|2008|) 



fearpov et all fl2008bh 



^Telinek et ajj fl2008h 



Hentunen et al.1 f|2008l ) 



J Novakl (|2008|) 



Kruglv et al.1 (120081 ) 



Tanvir et all fl2008h 



Levan et al.1 fj2008h 



n Converted from F606W. 
°Converted from F804W. 

v Calibrated to y-band and plotted as V in our figures. 
x No GCN observations were used in modeling. 
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Table 9. Bright SDSS standard stars 



RA 


Dec 


B 


V 


R 


I 


deg 


deg 


mag 


mag 


mag 


mag 


217.915229 


36.31513 


15.940 ± 0.010 


15.293 ± 0.006 


14.922 ± 0.008 


14.551 ± 0.007 


217.908143 


36.29125 


17.716 ± 0.012 


16.350 ± 0.006 


15.523 ± 0.008 


14.821 ± 0.007 


217.943617 


36.31667 


17.993 ± 0.012 


16.834 ± 0.006 


16.162 ± 0.008 


15.590 ± 0.007 



Table 10. Faint SDSS standard stars 



RA 


Dec 


g 


r 


i 


z 


deg 


deg 


mag 


mag 


mag 


mag 



217.908143 


36.29125 


17. 


,091 


± 





.005 


15. 


,816 


± 





.004 


15.310 


± 


0.004 


15.067 


± 


0.005 


217.943617 


36.31667 


17. 


,425 


± 





.005 


16. 


,409 


± 





.004 


16.058 


± 


0.004 


15.872 


± 


0.007 


217.871828 


36.29380 


17. 


,364 


± 





.005 


16. 


,903 


± 





.005 


16.769 


± 


0.005 


16.689 


± 


0.010 


217.894280 


36.34851 


18. 


.298 


± 


0, 


.007 


16. 


.955 


± 


0, 


.005 


16.370 


± 


0.005 


16.036 


± 


0.007 


217.902832 


36.33332 


18. 


,084 


± 





.006 


17. 


,741 


± 





.006 


17.626 


± 


0.007 


17.565 


± 


0.017 


217.981086 


36.29605 


18. 


.223 


± 


0, 


.007 


17. 


.866 


± 


0, 


.006 


17.746 


± 


0.008 


17.723 


± 


0.041 


217.931339 


36.34208 


19. 


.812 


± 


0, 


.014 


18. 


.342 


± 


0, 


.008 


16.818 


± 


0.005 


16.016 


± 


0.007 


217.931216 


36.27044 


19. 


.172 


± 


0, 


.010 


18. 


.382 


± 


0, 


.008 


18.110 


± 


0.009 


17.988 


± 


0.021 


217.918619 


36.31995 


20. 


,742 


± 





.026 


19. 


,264 


± 





.013 


18.483 


± 


0.011 


18.055 


± 


0.024 


217.913905 


36.28612 


20. 


.736 


± 


0, 


.027 


19. 


.266 


± 


0, 


.013 


18.253 


± 


0.010 


17.687 


± 


0.019 


217.928915 


36.33948 


20. 


,850 


± 





.029 


19. 


,462 


± 





.017 


18.214 


± 


0.010 


17.578 


± 


0.018 


217.928801 


36.30674 


20. 


,844 


± 





.028 


19. 


,471 


± 





.014 


18.689 


± 


0.012 


18.315 


± 


0.029 


217.955190 


36.33084 


19. 


,708 


± 





.014 


19. 


.506 


± 





.015 


19.505 


± 


0.023 


19.315 


± 


0.065 


217.926000 


36.33788 


20. 


.047 


± 


0, 


.016 


19. 


,640 


± 


0, 


.016 


19.474 


± 


0.022 


19.437 


± 


0.071 


217.932586 


36.27256 


21. 


.248 


± 


0, 


.033 


19. 


,710 


± 


0, 


.016 


18.342 


± 


0.010 


17.551 


± 


0.016 


217.946853 


36.30972 


21. 


.362 


± 


0, 


.040 


19. 


,909 


± 


0, 


.020 


19.003 


± 


0.016 


18.439 


± 


0.033 


217.971356 


36.30059 


21. 


.318 


± 


0, 


.039 


19.937 


± 


0, 


.020 


19.043 


± 


0.016 


18.571 


± 


0.035 



